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ABSTRACT 

We apply the reconstruction technique of Nusser & Dckcl (1992) to the recently 
available Point Source Catalogue Redshift Survey (PSCz) in order to subtract the 
phase correlations that are expected to develop in the mild non-linear regime of grav- 
itational evolution. We study the evolution of isodensity contours defined using an 
adaptive smoothing algorithm in order to minimize the problems derived from the 
non-comutivity of operators. We study the genus curves of these isodensity contours 
and concentrate on the evolution of the amplitude drops, a meta-statistic able to quat- 
ify the level of phase-correlation present in the field. In order to test the method and to 
quatify the level of statistical uncertainty, we apply the method to a set of mock PSCz 
catalogues derived from the N-body simulations of two 'standard' CDM models, kindly 
granted to us by the Virgo consortium. We find the method to be reliable in recover- 
ing the right amplitude drops. When applied to PSCz the level of phase correlations 
observed is very low on all scales ranging from S/i^Mpc to 60/i -1 Mpc, providing 
support to the theory that structure originated from gaussian initial conditions. 
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INTRODUCTION 



The IRAS Point Source Catalogue Redshift Survey ( Saun ■ 
ders et al. 199S| ) lends itself to topological studies because 



of its high sampling densities and large volume. It countains 
approximately 15000 galaxies to the full depth of the PSC 
(O.GJy). Its sky coverage is 84.1 per cent, where only the 
zone of avoidance is excluded, here defined as an infrared 
background exceeding 25MJysr _1 at lOO^tm, and a few un- 
observed or contaminated patches at higher latitude. The 
excluded regions are coded i n an angular mask, as shown 
Canavezes et al. (1996) analysed the topology 



Fig. 



of PSCz and showed that it is consistent with the topol- 
ogy of CDM models: The genus curves retain the w-shape 
characteristic of random-phase density fields even at small 
smoothing lengths, where non-linear evolution has already 
generated significant skewness of the 1-point PDF. These 
non-linearities are, however, detected by a depressed am- 
plitude of the genus curves - amplitude drops - consistent 
with those detected for the CDM models. Above 10 /i _1 Mpc 
strong phase correlations were not detected in PSCz and 
those that were detected at smaller scales are expected in the 
framework of mildly non-linear gravitational evolution. This 



supports the hypothesis that structure grew from random- 
phase initial conditions. 

In this paper we attempt to subtract these phase corre- 
lations by reversing gravity. If it is true that structure does 
indeed originate in random-phase fluctuations, then revers- 
ing gravitational evolution should provide us with an initial 
density fluctuation field free from phase correlations on any 
scale. 

Ther e are several methods to achieve this goal. Nusser 
& Dekel (1992), showed how to express the Zel'dovich ap- 
proximation in a set of Eulerian coordinates and reverse it to 
obtain a reconstructed density field. They tested the method 
with N-body simulations and found satisfactory results. Re- 
cently, some improvements have been made on the origi- 
nal m ethod (e.g. Gramann (1993)). Nusser , Dekel & Yahil 
( 1995| ) applied the Nusser & Dekel (1992) approximation 
to the 1.2 Jy IRAS redshift survey to recover the 1-point 
probability distribution function of the initial density field. 
Their results were consistent with Gaussian initial condi- 
tions. The PSCz presents us with an unprecedented number 
of resolution elements and should constitute therefore the 
best available data set on which to apply a reconstruction 
technique. 

The method we choose to follow here is the quasi-linear 
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Figure 1. Sky coverage of the PSCz survey in an aitoff projection. The dots represent the galaxies in the survey, the shaded regions are 
unobserved and comprise the angular mask. The galactic centre lies in the centre of the plot. 



method proposed by Nusser & Dekel ( |L992 ). In section [l] 
we describe this method; In section || we describe the con- 
struction of the density maps used as the input of the time- 
machine and describe our error estimates; Section ^ is de- 
voted to testing the method using N-body simulations of a 
CDM model and in section ^ we apply the method to the 
PSCz and present our results. 



1 THE RECONSTRUCTION METHOD 

The reconstruction method we follow is based on two 
premises: The velocity field is, when smoothed on a scale 
of a few Mpc, irrotational, and the Zel'dovich approxima- 
tion is accurate over the mildly non-linear regime. 

Our first step is to try to express the Zel'dovich approx- 
imation in a set of Eulerian coordinates. 

Given an initial comoving position q for a given particle, 
a final position x(q,t) will have the form 



x(q,t) =q + P(q,t). 



(1) 



The Zel'dovich approximation states that the displacement 
term P(q, t) can be written as a product of two functions, 
each a function only of one of the variables q or t, i.e., we 
can separate the variables q and t: 



x(q,t) = q + D(t)¥(q) 



(2) 



This is simply a linear approximation with respect to the 
particle displacements rather than density. Particles in the 
Zel'dovich approximation follow straight lines: 

<7x 



v(q,t) = o(i)— = o(t)£>(t)*(q). 
The density fluctuation is given by 



(3) 



6(q,t)=p/J(q,t)-l, 



(4) 



where J(q, t) is the Jacobian of the coordinate transforma- 
tion q — > x. 

As long as second order terms in S and v are negligible, 
S(t) oc D(t), where D(t) is the growing m ode so lution in 
linear perturbation theory (see e.g. Peebles (1980)). 

Unlike linear theory, where the density at a given posi- 
tion evolves according to the linear growth rate, under the 
Zel'dovich approximation infinite density can develop in a 
finite time as a result of the convergence of particle trajec- 
tories into a pancake, in order to in tegrate a density field 
back in time, Nusser & Dekel (1992) found the differential 
equation in Eulerian space which contains the Zel'dovich 
approximation. This is derived from the standard equation 
of motio n of d ust particles in an expanding Universe (e.g. 
Peebles (|l980|)): 



-tt + Hw = V$ s 

dt a 



(•») 



Here v stands for the peculiar velocity of the particle, i.e., 
v = ra(t)4r, where x are the comoving positions; V = 
(d/dx,d/dy,d/dz) and $ 9 is the gravitational potential 
which is related to the local density fluctuation via the Pois- 
son equation 



V 2 $ a = ^H 2 Q,a 2 5. 



(6) 



0(x,t) 

<t>g{x,t) 



Using the normalized variables 9 and tp g defined thus 
v(x,t) 



aD 
$g(*,t) 
a 2 _D 



(7) 
(8) 
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equations |B| and ^| reduce to 
3HQ 



> + 



2/(Q) 



(9) 



where f(Q) = D/HD ~ Q. os ( Peebles 1980[ ). 

Notice that the Zel'dovich approximation is contained 
in which is to say 

9 = (10) 

along the trajectory of a given particle, i.e., along a line 
of constant q. Hence, ^| is reduced, under the Zel'dovich 
approximation, to 

3HQ , 



2/(0)' ~ V ^- 
is therefore the gradient of a potential: 

= -V<p v , 
with 

, _ 2/(fi) 



3iTO 



<p„ + F(jt). 



(11) 



(12) 



(13) 



We can obviously set F(t) = since is the only physical 
(measurable) quantity. 

Equation HQ can now be expanded in Eulerian coordi- 
nates: 



d9 n 

— = O 

df 



at dt 



?(x,t) = 



O — + D(t)^(q) • Vfl = 
& ^+Z>(0- V)0 = O. 



(14) 



Making use of the identity (0 • V) = 1/2 (V0) 2 -0A(V A 0), 
taking into account the irrotationality of 0, and substituting 
|l2] we arrive at 

dip v , D 



V^ + ^(V0f 







dt 



dip v b 



ir+2 (v ^ 



Again, because only gradients of ip„ have a physical meaning, 
we can take F(t) to be zero. We finally arrive at the equation: 



(16) 



This is the differential equation which expresses the 
Zel'dovich approximation in Eulerian coordinates. Knowing 
<p v at any time enables us to compute ip v at any other time 
simply by integrating [l(] forwards or backwards. 

The first step is to calculate tp v at the present time. 
That can be easily achieved using Poisson's equation, which 
in Fourier space reads: 

- k 2 $ 3 = ^H 2 na 2 ~S, (17) 

where the tildes denote the Fourier transforms. Using [] we 
arrive at: 



Dk 2 <p v = 8. 



(18) 



So, in order to obtain the velocity potential of a given 
smooth density field with periodic boundary conditions, we 
first transform it to Fourier space using some FFT code, 



Table X. Parameters of the selection function of PSCz. 



n qqi +0.068 
u,ro -0.073 



P 



3.445 



+0.173 



7 



1.925 



-0.162 
-0.153 



z* tp [^Mpc" 3 ] 

0.025341°;°°^° (141.3 ± 2.4) x 10" 



we then divide the obtained field by — k 2 (normalizing D to 
unity at present), and then transform back to real space to 
obtain <p v . <p v is then integrated back in time in the most 
trivial way: At every step, we calculate 1/2 (V<p v ) 2 from the 
potential field; this is then used as a first order Taylor cor- 
rection to predict the value of <p v at the next step. Once 
a suitable number of integration steps has been performed, 
we use equation [l8] again to obtain the value of 8 at the 
corresponding value of D. 

For our application - topology, we expect the exact num- 
ber of integrations to be irrelevant. After a certain number of 
iterations we expect to reach the linear regime in all scales, 
which means that the genus curves of the density fields will 
no longer change. This particular characteristic of the genus 
can in fact be useful to test the convergence of the method. 
We obtained several reconstructed density fields for differ- 
ent values of the number of integration steps, and found the 
topologies to converge. This result is shown is Fig. ^. 



2 CONSTRUCTION OF THE DENSITY MAPS 

The construction of the density maps we follow is, to a large 
exte nt, e quivalent to the method followed by Canavezes et 
al. (199c). We employ the same fit for the PSCz selection 
function s(z): 



8{Z) 



(i + G 



\-r\P/-y' 



(19) 



with the parameters shown in Table jSpringel fe Whi te 
199c), and estimate the density p(r) by 



p(r) tx 



m(r) 
s{r) ' J 



(20) 



where m(r) is the discrete point distribution. However, 
the angular mask needs to be treated with care. Because the 
structure we can see behind the mask will have a significant 
effect on the evolution of structure on the whole observable 
area, we need to resort to some sort of filling, prior to using 
the time-machine. 

We consider two different types of filling: One, a random 
filling, in which fake objects are placed randomly over the 
masked region with an average number density equal to the 
average number density of the whole survey (weighted by 
the selection function); And another, a cloning, in which 
fake objects are placed randomly in each bin, but with an 
average number density equal to the average number density 
of the neighbouring observable bins, again weighted by the 
selection function. 

On both cases, we create a box containing a sphere 
of radius R ma x as defined in Table ^| R m ax is the maxi- 
mum distance up to which the average distance between two 
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Table 2. The smoothing lengths adopted for the topological anal- 
ysis of the PSCz survey. Listed are the adopted survey depth 
fimai, the resulting number N rcB of resolution elements and the 
number -/V gal of galaxies inside the survey volume. 



ivipcj 


-Kmax [ H IVlpCJ 


AT 

JVres 


AT 


F 






zlo.o 




a 



4 1 . ID 


QfV7 Q 


OOOU 


7 


58.35 


366.3 


4928 


8 


68.58 


398.5 


5909 


10 


86.91 


415.3 


7510 


12 


103.19 


402.3 


8775 


14 


118.04 


379.2 


9681 


16 


131.81 


353.7 


10356 


20 


156.96 


305.8 


11309 


24 


179.79 


266.0 


11968 


28 


200.97 


233.9 


12496 


32 


221.79 


210.6 


12833 


40 


263.44 


180.7 


13339 


48 


305.09 


162.4 


13684 


50 


315.50 


158.9 


13748 


56 


346.73 


150.2 


13924 



neighbouring galaxies in PSCz is smaller than the adopted 
smoothing length. The size of the box will then be -j=Rmax. 
Although boundary conditions are not periodic, we smooth 
the density field in Fourier space assuming periodic bound- 
ary conditions. We consider this to be preferable to zero 
padding, as this would create an artificial boundary that 
would eventually affect the gravitational time-machine. 

We then apply the Zel'dovich time-machine on the box 
thus obtained for both filling techniques. In the subsequent 
topological analysis, we limit ourselves to the sphere inside 
the box with radius R ma x ■ 

Smoothing Procedure 

Ideally, one would want to apply the operator time-reverse 
directly upon the unsmoothed density field obtained from 
PSCz, and only then smooth it on a range of scales to 
be able to calculate its genus curve. However, in order to 
obtain meaningful results, we need a smooth field a pri- 
ori. This poses a very important problem: These operators 
(time-reverse and smoothing) do not in general commute. 
In other words, the topology of the final density field can 
be significantly different whether the smoothing operation 
is performed before or after applying the time-machine to 
our original field. 

To illustrate this point, let us consider a Gaussian ran- 
dom density field smoothed on some scale A. After applying 
the time-machine this density field will still be Gaussian. 
However, if we are to naively calculate its genus curve choos- 
ing for isodensity contours the same isodensity contours de- 
fined prior to applying te time-machine, we would wrongly 
conclude that the field is not Gaussian. In order to circum- 
vent this problem we would be forced to smooth the final 
field once again on a larger scale, thus reducing considerably 
the statistical significance of our results. 

There is, however, an alternative way to solve this 
dilemma: One can try to find a smoothing operator that 
comutes with our reconstruction operator. 



How can we look for such a smoothing operator? One 
way of ensuring this is to find an operator that does not 
change the density field if applied once again at a latter 
stage. By other words, if we ensure that our final density 
field remains the same regardless of how many times we per- 
form the smoothing operation throughout the time-reversing 
process, then it is because these two operators do indeed co- 
mute. A class of such operators are the adaptive smoothing 
operators. By adaptive smoothing, we understand a local 
smoothing of variable smoothing length according to the lo- 
cal structure. There are several types of adaptive smoothing. 
The ideal adaptive smoothing operator ensures that the to- 
tal mass enclosed by isodensity contours remains constant 
throughout the time-reversing process. A simplified version 
assumes spherical symmetry. At each point we employ a 
smoothing length A such that A 3 cx 1/p = 1/ (p + S). The 
proportionality constant defines the characteristic smooth- 
ing length, which is the smoothing length employed when the 
overdensity 8 vanishes. For highly dense regions, a smaller 
value of A will be chosen, whereas voids will be smoothed 
with a larger value of A. 

When we apply the time-machine operator on a map 
smoothed in this way, the isodensity contours around a clus- 
ter will move, but the total mass enclosed by them will re- 
main approximately constant, as long as the spherical sym- 
metry hypothesis is a good approximation. 

In our application to the PSCz we use a spherically sym- 
metric adaptive smoothing algorithm. We start by creating 
a set of maps obtained from the original map (this might 
be either the PSCz data or simulation data) by smoothing 
it on a range of scales around a given characteristic scale l , 
using Gaussian filters of the form 



where A varies around lo- We then look for the appropriate 
value of A at each position x by enforcing the equation 

(l + 5 x )\ s = ll (22) 

where 8\ is the density contrast when the original map is 
smoothed on a scale A. This ensures that the mass enclosed 
on a sphere of radius A is, to first order in a Taylor expansion, 
constant throughout the whole map. 

Error Estimates 

There are three different sources of error that enter our re- 
sults: Shot noise, cosmic variance and the errors associated 
with the Zel'dovich approximation itself. 

The most accurate and realistic way of estimating these 
errors, which ultimately enter the genus curves of the re- 
constructed PSCz density fields, is to analyse the topology 
of fake PSCz catalogues drawn from N-body simulations of 
some standard CDM model. Since we intend to test the valid- 
ity of the Zel'dovich time-machine with N-body simulations, 
it seems reasonable to extent this philosophy to the calcula- 
tion of the statistical errors themselves. This is achieved by 
seeking a galaxy number density field with a Poisson distri- 
bution, whose expectation value is identical to the density 
field of the N-body simulation multiplied by the PSCz se- 
lection function. This is equivalent to observing the density 
field of the N-body simulation in a similar way to the way 
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Table 3. Parameters of the examined CDM models. The simula- 
tions have been done by the Virgo collaboration. 





rCDM 


ACDM 


Number of particles 


256 3 


256 3 


Box size[ /i _1 Mpc] 


239.5 


239.5 


Zstart 


50 


30 


Qo 


1.0 


0.3 


n A 


0.0 


0.7 


Hubble constant h 


0.5 


0.7 


r 


0.21 


0.21 


eg 


0.60 


0.90 


Mass per particle [1O 1o /i _1 M0] 


22.7 


5.8 



PSCz observes the density field of the real Universe. This 
galaxy number density field is then divided by the PSCz se- 
lection function again to obtain a "PSCZ-noisy" distance- 
independent estimate of the real galaxy number density. 
Each of the mock PSCz catalogues is adaptively smoothed 
using the algorithm described above and then used as input 
in the Zel'dovich time-machine. The genus curves of the re- 
constructed fields are calculated and the variance obtained 
over 10 mock PSCz catalogues is used as the statistical error 
estimate. 



3 TESTING THE METHOD WITH N-BODY 
SIMULATIONS OF CDM MODELS. 

As we mentioned previously, we intend to test the regime 
of validity of the Zel'dovich time-machine by means of an 
N-body simulation of a CDM model. We apply the time- 
machine on the present density field drawn from the sim- 
ulation and compare the density field thus obtained with 
the density field drawn from the original test-mass positions 
used as input in the simulation. More specifically, we need 
only to compare the topologies of the reconstructed and orig- 
inal fields through their genus curves. 

For this purpose, we use two N-body simulations corre- 
sponding to two cold dar k matter models, k i ndly provided 
by th e Virgo consortium ( |Colberg et al. 1997]; |Jenkins et al 



1997). The simulations hav e been performed with an A P^M 
SPH code named HYDRA (|Couchman fc Pearce 1995[) . Here 



we consider the SCDM model and the rCDM model, which 
parameters are shown in Table ^. Because this simulations 
contain CDM on periodic boxes of size 239.5 /i _1 Mpc and 
use such a large number of particles, they constitute an ideal 
ground for this test. 

We start by binning the particles in cells of side 
length 2 ft _1 Mpc and smooth further on some characteristic 
smoothing length A, using the adaptive smoothing algorithm 
described in section ^, in an analogous procedure to the way 
we treat the PSCz galaxies. This "double" smoothing is re- 
quired in the real Universe in order to minimize shot noise 
effects, but it is also necessary in order to eliminate severe 
non-linearities since the Zel'dovich time-machine is expected 
not to work over such regimes. It is also needed to smooth 
out regions of orbit-crossing. 

The next step is to calculate the genus curves of both 
the original density field that is used as input in the simu- 
lation, smoothed in the same manner as described above, 
and the reconstructed density field obtained after apply- 



ing the Zel'dovich time-machine to the present density field 
smoothed in the same way. In order to compare both genus 
curves we compute their amplitudes and their amplitude 
drops. It is very important to note that if we are to ob- 
tain gaussianized versions of the density fields at the present 
time, they will not have, in principle, the same genus am- 
plitudes as the density fields themselves, even when these 
are in the Gaussian regime. This is so because of the very 
nature of adaptive smoothing. When a given density field 
is smoothed adaptively its genus curve will have a different 
shape and amplitude than the shape and amplitude of the 
genus curve of the same field when this is smoothed with a 
constant Gaussian window. We need to be extremely careful 
not to draw the wrong conclusions about the Gaussian or 
non-Gaussian nature of our density fields. After applying the 
time- machine, however, the amplitude of fluctuations will be 
reduced significantly. In fact they will be reduced to such an 
extent as to make the field look almost homogeneous. This 
means that the isodensity contours on adaptively smoothed 
maps will be very close to the isodensity contours on maps 
smoothed using a constant Gaussian window. Hence, it is 
possible to calculate amplitude drops and determine the 
Gaussian (or non-Gaussian) nature of our density maps af- 
ter applying the Zel'dovich time-machine, i.e., to determine 
the Gaussian (or non-Gaussian) nature of the reconstructed 
fields. 

In order to test the convergence of the reconstruction 
method, we obtained several reconstructed density fields for 
diffent values of the total number of integration steps, as 
mentioned in |l|. Fig. |^ shows a particular slice of the recon- 
structed density fields at different integration steps when the 
method is applied to the maps at z = obtained from the 
rCDM simulation by smoothing on characteristic lengths of 
8/i _1 Mpc and 10/i _1 Mpc. As it is evident from this Fig- 
ure, the topology of the density fields is undistinguishable 
for values of z greater than 4 on smoothing scales of both 
8 h~ 1 Mpc and 10 ft _1 Mpc. Only the amplitude of the density 
fluctuations changes, indicating that we are now in the linear 
regime. In Fig. ^ this is made even clearer. Here we show a 
point-by-point comparison of the reconstructed rCDM den- 
sity fields at redshifts z = 4 and z = 9, for the characteristic 
smoothing lengths of 8/i _1 Mpc and 10/t -1 Mpc. We only 
show one in eight of all points, chosen randomly. It is ob- 
vious from this plot that the shape of fluctuations did not 
change from z — 4 to z = 9, for either of the smoothing 
lengths adopted. It is also obvious that the amplitude of 
fluctuations was reduced by a factor of 2, as it is expected 
in the linear regime (notice that in the linear regime the 
growing mode of fluctuations varies as 1/(1 + z)). Hence, 
our time-machine shows the correct assimptotic behaviour. 



Results 

Figure ^| shows the isodensity contours obtained for the 
SCDM simulation and the rCDM simulation, when the 
fields are adaptively smoothed on a characteristic scale of 
5/i _1 Mpc in a comoving box of size 240/i _1 Mpc. The first 
row shows the fields at the present time. The second row 
shows the fields after applying the Zel'dovich time-machine 
back to a redshift of 10 and the third row shows the orig- 
inal fields used as input in the simulations. In each panel 
the thick line represents the contour where the density con- 
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Figure 2. A particular slice of the reconstructed fields for the tCDM simulation for different values of z. The left panels show the fields 
smoothed adaptively with a characteristic smoothing length of 8/i _1 Mpc, whereas the panels on the right show the fields smoothed with 
a characteristic smoothing length of fO/i _1 Mpc 



tCDM 

0.15 - A = 8/T'Mpc 




-o m i i . 

-0.20 -0.10 0.00 



tCDM 

X = lOTT'Mpc 



-0.10^ 
0.40 -0.2 




Figure 3. Point-by-point comparison of the reconstructed tCDM density field at redshift z = 4 and the reconstructed rCDM density 
field at redshift 2 = 9, when the adopted characteristic smoothing lengths are of 8/i -1 Mpc and f0/i -1 Mpc. 
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TCDM_origina 



,50 L/> 



U -- 



\ t 



Figure 4. Isodensity contours of the tCDM simulation. The characteristic smoothing length adopted was 8 h ^^Mpc. The separation 
between contour levels is 0.2/(1 + z). The zero contour is heavy; positive contours are solid; negative contours are dashed. 



trast is zero. Dashed lines represent contours of negative 
density contrast whereas solid lines represent positive den- 
sity contrasts. The contour spacing is 0.2/(1 + z). Because 
this is normalized to the linear reconstruction case where 
S oc 1/(1 + z), the linear theory reconstruction density maps 
look exactly like the present day maps (first row), albeit 
their much lower fluctuation amplitude. From the contours 
on both the reconstructed maps and the original maps we 
notice that their fluctuation amplitudes are larger than the 
fluctuation amplitude of a linearly reconstructed map. This 
result is expected because fluctuations will grow faster at the 
later stages of evolution, in the mildly non-linear regime. It 
is obvious from Fig. ^ that the Zel'dovich time-machine is 
able to change the rank order of isodensity contours. This 
means that the genus curves will also be changed. In fact we 
can predict, just by looking at Fig. ^ that the amplitude of 
the genus curves will increase after applying the Zel'dovich 
time-machine, as it is expected. 

Figures ^ and ^ show the genus curves obtained for the 
SCDM model and the rCDM model respectively, at selected 
characteristic smoothing lengths, namely at 5/i~ 1 Mpc, 
8ft _1 Mpc and 14/i -1 Mpc. We restrict our analysis to this 
range of smoothing lengths, as we know a priori that the 
Zel'dovich time-machine is only relevant in the mildly non- 
linear regime. In both figures, the first column shows the 
genus curves of the reconstructed field and its randomized 
counterpart, whereas the second column shows the genus 



SCDM_z=0 
tCDM_z=0 
SCDW_recom 
TCDM_recon: 



Figure 7. Amplitude drops measured for the reconstructed 
SCDM simulation and rCDM simulation. Also show for com- 
parison are the amplitude drops measured for the simulations at 

2 = 0. 



curves of the original fields, used as input for the simula- 
tions. We deliberately do not show the genus curves of the 
present day adoptively smoothed density fields, since, as we 
argued in the previous chapter, the isodensity contours are, 
at this stage, different from the isodensity contours of the 
fields smoothed using a constant kernel. 
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SCDM_reconstructed X=^/T'Mpc SCDy_oriqina A=5/r'Mpc 






SCDM_reconstructed X=14/r'Mpc 5CDM_oricingl A=14/T'Wpc 




Figure 5. Genus curves for the SCDM model, at selected smoothing lengths. The left column shows the genus curves of the reconstructed 
field (thick solid lines) and of the randomized version of that field (thin solid lines), whereas the right column shows the genus curves of 
the original density field used as input for the N-body simulations. The dashed and dotted lines are the best fit random-phase curves to 
each of the fields. 



It is striking in Fig. g and Fig. g the consistency be- 
tween the genus curves of the reconstructed field and of its 
randomized counterpart. This is particularly true for the 
genus amplitudes, which indicates that there is very little 
phase correlation in the reconstructed fields on all scales, as 
it is expected from Gaussian initial conditions. With regard 
to the genus amplitudes per se, we find consistency with 
the genus amplitudes of the original density fields, although 
we notice a tendency for these to be slightly higher in the 
reconstructed fields. These means that the second moment 
of the power spectrum is recovered slightly in excess. This 
becomes more important as we approach the smaller scales, 
which is to say as the system becomes more non-linear. 

It is also interesting to notice a depression at high val- 
ues of v on the genus curves of the reconstructed fields. We 
find this depression to be present in all cases, both for the 
Virgo simulations and the PSCz data. We expect this to be 
a particular feature of the reconstruction technique, i.e., ul- 
timately dependent on the Zel'dovich approximation itself. 
At high values of v the regime is considerably non-linear 
and so we expect the Zel'dovich approximation to perform 
rather poorly. This problem is enhanced by the fact that the 



density fields have been smoothed adaptively. Nevertheless, 
we do not expect this feature to affect our results consider- 
ably as the calculation of the amplitude of the genus curves 
is restricted to the range —\<v<\. 

Note that there is no sampling noise in these density 
fields. The tremble in the genus curves is ultimately due to 
cosmic variance. 

In Fig. ^ we plot the amplitude drops obtained for 
the reconstructed SCDM model and for the reconstructed 
rCDM model as a function of smoothing length. Also shown 
are the amplitude drop s obtained for the models at z = 
( |Canavezes et al. 1998 ). In all scales, the values obtained 



for the reconstructed fields are closer to unity than those 
of the present-day fields. This indicates that the Zel'dovich 
time-machine is a good tool in recovering Gaussian initial 
density fields. However, at small scales, we detect a slight 
departure from the expected value of unity, even for the re- 
constructed fields. This is not surprising since at this scales 
non-linearities become strong. The degree of strength of the 
non-linearities depends on the particular model. As we see 
from Fig. [7] they are more important in the rCDM model 
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Figure 6. Genus curves for the tCDM model, at selected smoothing lengths. The left column shows the genus curves of the reconstructed 
field (thick solid lines) and of the randomized version of that field (thin solid lines), whereas the right column shows the genus curves of 
the original density field used as input for the N-body simulations. The dashed and dotted lines are the best fit random-phase curves to 
each of the fields. 



than in the SCDM model , in agreement to what was found 
by Canavezes et al. (1998). 

Fig. ^ shows the genus amplitudes of the reconstructed 
fields for both the SCDM model and the rCDM model, to- 
gether with the genus amplitudes of the original fields. In 
the case of the rCDM model the agreement between the 
original genus amplitude and the genus amplitude of the re- 
constructed field, is striking, in particular when we consider 
scales above ~ 7/i~ 1 Mpc. For the SCDM model the agree- 
ment is not as good. The reconstructed genus amplitudes 
appear to be slightly higher than the true original ones. 

These results indicate that the Zel'dovich time-machine 
is effective in recovering the right genus amplitude drops on 
scales larger than ~ 8/i~ 1 Mpc, although recovered ampli- 
tudes cannot be considered reliable in the strict sense of the 
word. 



RECOVERING THE INITIAL DENSITY 
FIELD FROM PSCZ 



As mentioned previously, one of the difficulties in using the 
PSCz in an attempt to recover the initial density fluctua- 
tions in the Universe, is the fact that the masked region 
can make up to ~ 20% of the whole observable area. Since 
this region will have a definite gravitational influence on the 
observable area when we attempt to reconstruct to original 
density field, some sort of filling is essential. As mentioned 
in ^ we employ two different techniques: The random filling 
and the cloning. Fig. |To| shows how the two different fillings 
appear to the eye. On the upper panel we show the map 
where the mask has been filled randomly and on the lower 
panel we show the map where the mask has been cloned. 
Although it is impossible to the eye to recognize any signif- 
icant difference, we will be able to detect some noticeable 
differences on the topologies of the reconstructed density 
fields. 
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Figure 9. Genus curves of the reconstructed PSCz fields. 



Results & Discussion 

Fig. |^ shows the genus curves obtained for the reconstructed 
PSCz fields. On the first column we plot the genus curves of 
the fields for which the mask has been filled using a cloning 
technique, and on the second column we plot the genus 
curves of the fields where the mask has been filled randomly. 
The thick solid lines refer to the reconstructed PSCz fields 
proper; whereas the thin solid lines refer to the randomized 
versions of those fields, i.e., the fields obtained from the re- 
constructed density fields by randomizing phases in Fourier 
space subject to the reality constraint <5k = <5_I k and keep- 
ing the same power spectrum. The dotted and dashed lines 



are the best fitting random-phase curves to both the recon- 
structed PSCz density fields and their randomized versions. 

It is striking to notice the proximity between genus 
curves, even at small smoothing lengths. The statistical er- 
rors drawn from the N-body simulation are not shown in 
Fig. ^ because of the high degree of correlation between the 
points in the genus curves. 

Independently of whether the mask has been randomly 
filled or cloned, the recovered genus curves seem consistent 
with random-phase Gaussian fluctuations because the am- 
plitude drops appear small. However, the amplitude them- 
selves seem to depend on the mask filling technique. 

In Fig. ^ we plot the amplitude drops obtained for the 



© 0000 RAS, MNRAS 000, 000-000 



Recovering the Topology of the Initial Density Fluctuations Using the IRAS Point Source Catalogue Redshift Surve 



PSCz random mask PSCz cloned mask 









...-£55..;,: r^&- . %~ "■ .. ... .is 




. . - V" i'-' . : ' * • ; '*.. '•*" E- ; ' ~ ' 

',35_., .„».*£;„ >+^",, ; . v .. .(}. . ', 4p -.CP :.r55E.:.l, , .18 

























Figure 10. AitofT projections of PSCz galaxies and filled masks. In the upper panel the mask is filled randomly whereas in the lower 
panel the mask is filled using a cloning technique. 
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Figure 11. The amplitude drop R measured for PSCz. Also shown are the results obtained by Vogeley et al.1994 for the CfA survey. 
Here the errors are taken to be the uncertainty Vogeley et al.1994 report for mock catalogues extracted from a LCDM model 
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Figure 12. The amplitude drops measured for the reconstructed PSCz density fields. In the upper panel the mask has been cloned, 
whereas in the lower panel the mask has been filled randomly. Also shown are the amplitude drops obtained for the reconstructed N-body 
simulations. 
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Figure 8. Genus amplitudes measured for the reconstructed 
SCDM simulation and rCDM simulation. Also shown are the 
genus amplitudes of the original fields. 



reconstructed PSCz density fields, as a function of smooth- 
ing length, as well as the error bars drawn from the rCDM 
model following the method outlined in chapter ^. The am- 
plitude drops of the reconstructed N-body simulations are 
also shown for comparison. 

In contrast with Fig. O, where phase correlations were 
found for PSCz at small smoothing lengths, the recon- 
structed PSCz fields do not show any significant phase corre- 
lations on scales ranging from 5 /i _1 Mpc to 14 /i _1 Mpc. This 
reinforces the hypothesis that density fluctuation originate 
from random-phase Gaussian initial conditions. 
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